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We discuss mean field theory of glasses without quenched disorder focusing on the justification of 
the replica approach to thermodynamics. We emphasize the assumptions implicit in this method 
and discuss how they can be verified. The formalism is applied to the long range Ising model with 
orthogonal coupling matrix. We find the one step replica-symmetry breaking solution and show that 
it is stable in the intermediate temperature range that includes the glass state but excludes very low 
temperatures. At very low temperatures this solution becomes unstable and this approach fails. 



The thermodynamics of glasses without quenched dis- 
order is a long standing problem in statistical physics. 
The interest to this problem was renewed recently when 
it was understood that powerful methods developed for 
the glasses with quenched disorder can be often applied 
to this problem. [0-^. In both systems the local magne- 
tization in the ground state varies from site to site and 
different sites are typically non-equivalent. The qualita- 
tive reason why glasses without quenched disorder are 
more difficult to describe theoretically than spin glasses 
is the following. The mean field theory has to operate 
with the average magnetization (or its moments), not 
with quantities which depend on a realization and a par- 
ticular state. The average quantities appear naturally in 
spin glasses after averaging over quenched disorder which 
makes all sites equivalent. 

A few methods were suggested to overcome this diffi- 
culty for the glasses without quenched disorder. First, a 
mapping of some glass models to the quenched disordered 
problems was suggested this method has an obvious 
disadvantage that such mapping is difficult to guess. Sec- 
ond, it was noted that a typical dynamics in a glassy 
system leads not to a ground state but to one of many 
metastable states providing an effective averaging mech- 
anism which makes all sites equivalent even for glasses 
without quenched disorder. This method has a disadvan- 
tage that dynamical equations are much more difficult to 
solve than the statical ones. Very recently the cloning 
method was proposed that is based on the idea that even 
at low T a system of m clones might be distributed in its 
phase space over many low lying metastable states if m is 
chosen correctly and the properties of all these states are 
essentially equivalent to those of the ground state. 
Generally, the partition sum of m weakly coupled clones 
is e~^'-™^^~'^'=°"^^^^\ where sum goes over free en- 
ergies (per site) of metastable states, F, and Sconf{F) 
is their configurational entropy {Sconf = ^^{■N' states))- 
Assuming that dSconf{F)/dF is finite at the lowest F as- 
sociated with the ground state one needs to chose m ocT 
at low T in order to avoid a complete dominance by a 
single (ground) state and the problems with site non- 
equivalence mentioned above. Distributing the system in 
the phase space provides the effective averaging mecha- 



nism in this approach. The main assumptions implicit 
in this approach are that low lying metastable states are 
not correlated (otherwise, averaging over them would not 
remove completely the non-equivalence of different sites) 
and that configurational entropy associated with these 
states behaves well as a function of energy at low ener- 
gies permitting the "right" choice of m. 

The goal of this paper is to provide an alternative the- 
oretical framework to the cloning method which, albeit 
somewhat similar in formalism, uses different physical ar- 
guments for its justification and allows to check the main 
assumption of the method mentioned above. We apply 
this method to the Ising version of the periodic long- 
range Josephson array model which is a simple example 
of a glass without quenched disorder and we show that in 
this model the main assumptions of the method are cor- 
rect in the intermediate temperature range but become 
wrong at very low temperatures. The main idea of the 
approach is that in a system with many low-lying states 
even a small random field is able to change the energy 
balance between the states and pull down a different state 
making it a new ground state of the system. Averaging 
over this random field is equivalent to the averaging over 
low lying metastable states. Specifically, in a spin system 
we add to the physical Hamiltonian a magnetic field part: 
Ti Ti -\-^^hiSi with small random hi. The resulting 
change in the energy of a typical metastable state is of 
the order oi^/Nh; because this energy interval contains 
a large amount of metastable states, we expect that a 
small non-zero field would result in a large rearrangement 
of their energies but would not change the properties of 
individual states. Averaging over the random field con- 
figurations is performed in the usual way introducing n 
replicas of the system and taking the limit n — > 0. The 
assumption of uncorrelated states is equivalent to one 
step replica symmetry breaking (IRSB) formalism; fur- 
ther, in this case this method is formally equivalent to 
cloned liquid approach if the size of the blocks in IRSB 
is equal to the number of clones (see |^ for the discus- 
sion of replica method vs clones for quenched disordered 
glasses). From the above discussion it is evident that 
another assumption implicit in this approach is that the 
energy spacing between low lying states should be much 
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less than 0{yN) if it is too big a small magnetic field 
will not be sufficient to rearrange low lying states, if it is 
too small, e.g. dSconf /dF\pg — oo, the effect of random 
field will be too large and no sensible limit hi ^ is 
possible. The latter situation seems to happen in the pe- 
riodic long range Josephson array with flux 2it per strip 
1^ when all states are exactly degenerate and Sconf{F) 
is very singular at T = 0. 

We now provide the details of our formalism and its ap- 
plication to the simplest model of a glass without disor- 
der. Our model consists of two sets of Ising spins (which 
we shall refer as "upper" and "lower" in the following) 
interacting via 
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(1) 



Here the spin Sim has a site index (m = 1 . . . N) and a 
components index i = 1,2 corresponding to the upper 
and lower spins and matrix J is 



J 777,7), 



J rji7 

J mn 



(2) 



with J„„ - (Jov^7]V)cos(^(m- f/2)(n- f/2)). For 
a = 1/2 we obtain the orthogonal limit "Y^^JmnJnk — 
J^Smk, in what follows we shall focus on this case. This 
Ising spin model is similar to the XY spin model of long 
ranged Josephson array and to the Bernasconi model 
1^. As well as in these models its lowest states corre- 
spond to "pseudorandom" sequences with flat Fourier 
transform. So, we expect that this model also displays 
glassy properties, in particular that it has extensive con- 
figurational entropy at low temperatures. Further, one 
expects that in a model with long range interaction the 
barriers separating metastable states become infinite in 
the thermodynamic limit. We have verified numerically 
that the configurational entropy in this model is indeed 
extensive and its dependence on energy is similar to the 
one obtained for other infinite range glasses (see Fig. 
1). Note, however, the important difference between this 
model and the XY spin model of [||: in the orthogonal 
limit the ground state of Ising model does not become 
extensively degenerate (i.e. degeneracy stays finite as 
TV — > oo, see Fig. 1) whereas in the XY spin model 
the ground state becomes extensively degenerate in the 
unitary limit making it very complicated 

Taking the Gaussian distribution for the random mag- 
netic field {hi hj) 
nian 



2 /iq 6i,j we get the replica Hamilto- 



a.p.i 



(3) 



where the replica indexes a, (3 run from 1 to n and the 
limit n — > should be taken. The glass transition cor- 
responds to an appearance of a non replica-symmetric 
solution of the Hamiltonian (|^) in the limit ho 0. 
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FIG. 1. Main panel: configura- 

tional entropy Sconf = \n{J^)/N at T = as a function of 
state energy obtained from direct numerics on systems up to 
N = 27 size. Insert: total Sconf and the degeneracy, Ng, of 
the lowest energy state at T = as a function of system size. 

In the large N limit a long range model containing N 
sites can be reduced to an effective single-site model with 
a free energy density 



(4) 



where S" is Ising spin field retaining only replica and 
component index dependence, B is an order parameter 
matrix in the replica space. The function ^{B) can be 
determined from the condition that all single site correla- 
tion functions of the model (^ coincide with the correla- 
tion functions of the original model (||). Instead of com- 
paring the spin correlation functions of these two models 
it is easier to decouple Ising spins by auxiliary field -0, 
sum over Ising spins and compare the correlation func- 
tions of conjugate field ip in the two new models 



-E^(^")- (6) 



where V{ip) = In 2 cosh(?/;). For both models one can 
construct a formal perturbation theory in the interaction 
ln2cosh(?/;Q -') and verify that these expansions coincide. 
We begin with the model (||) . Inspecting the terms of the 
perturbation theory for the correlator G"^ ^„ = {i^fm'^jn) 
one verifies that in the leading order in it is given 
by G = [TJ~^ — E]~^ with the self energy S which is di- 
agonal in the site index: S = A6mn Sij . This approach 
is similar to a locator expansion but in our case the 
locator A might be non-trivial in the replica space. Using 
the orthogonality of J we get that the single site corre- 
lation function Qap = G^f j„j (that we need to establish 
the correspondence between the models) becomes 
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(7) where X = Jq — + rjm)'^ and the function / is 



where jo = (3 J. 

Now we turn to the model (|^). Here the self energy 
is diagonal in the site index by construction, further, the 
interaction part of this model is the same as for model 
assuming that their single site correlation functions 
coincide we conclude that their single site sclf-encrgies 
are equal as well. Thus, the spin correlator obtained for 
this model isG = [B~^ — A]~^ , comparing this expression 
with we conclude that B — j^A . 

The correlator of the dual field ip can be related to 
the correlator of original spins: consider a Gaussian 
transformation leading to the model (||) exp{SBS/2) = 
J dip exp{—il;B^-^i/j /2 + Si/j) and use it to express G via 
correlator D^p = {S"'S^), we get: Q =B + BDB. Solv- 
ing this equation for the spin correlator D and using (j^) 
and the relation B = j^A we obtain 

D = B[f,-BT'- (8) 

Finally, the saddle point condition for the free energy (Q) 
2D = — 7'(i?), therefore integrating Eq.dq) we find 



7(S)-ln(l-Jo-2i?2), 



(9) 



Note that the free energy (Q) coincides with the free en- 
ergy of the model considered in Ref. although their 
properties at finite N are markedly different. Further- 
more, this free energy is the same as obtained by fidu- 
ciary Hamiltonian approach 

Paramagnetic state. In this state we take the replica 
symmetric ansatz i?Q,/3 = A*^c(,/3 ^-nd free energy (Q) be- 
comes 

T/T - [Info - HJo " A'')]/2 - /i - 2 In 2. (10) 



Variation with respect to /i gives /i = [\/T-l-4jf — l]/2. 
Usual thermodynamic relations between energy and en- 
tropy give E = -Tfi, S = In [4^/ jo]. 

One can see that the entropy of the normal solution 
becomes negative at T < Tk = Jo/(4%/T5) « 0.064550 Jo 
which is the Kauzmann temperature for this model ; 
one expects that the glass transition takes place at some 
temperature, Tc, above Tk- 

Glass state. At the glass transition temperature Tc 
the replica symmetry is broken, we assume that it is de- 
scribed by one step replica symmetry breaking (IRSB) 
and then verify that it is indeed a stable solution below 
Tc. The IRSB ansatz is i3a,/3 = M Sa^p+r] Ra^f3, where the 
matrix i? is a block-diagonal matrix consisting of m x m 
blocks with all elements equal 1, we get the free energy 
functional 



PJ^ = [log j2 _ (1 _ i/„,) In (j2 _ f,^)]/2 



21n2- (lnX)/2m-/i-2/(?7,m)/TO, (11) 



f{r],m) = In 



Pjn{z)dz 



P„,(z) = l^cosh"(zV^). 



/27r 



(12) 



Taking the derivatives of T with respect to /i, ry, m we get 

(13) 



m / Jo ^ M Xm 



(rym + /i)/X + q{m - 1) + 1 = 0, 



(14) 



2m? 



log [{jI-p?)/X] 



2_ _d_ 

m dm 



f{r),m) 



— r] [fi + rim,)/mX — 2 f{ri,m)/m'^ = 0, (15) 

where q — J tanh^(z)Pm(z)(iz/ J Prn{z)dz is the spin 
overlap of different replicas belonging to the same block 
^Q,/3 = i^^Q)Sa.f3 + Q Ra,0i that coiucidcs with Edwards- 
Anderson (EA) order parameter, 
solved with respect to to, ?y giving 



Eqs.(y,hJ) can be 



^ = 77 



1 + (1 — q + qm)rjm 

q/{^ - q) - 2v - q + qm) 



(16) 



and jg = + /i/(l — q). For a given to we can solve 
Eq. (n5|) numerically with respect to rj and get all quanti- 
ties as functions of to. The resulting dependence of to(T) 
for Jo = 1 is shown on Fig. 2. In the limit n ^ 0, 
the values of m should lie within the interval (0, 1) and 
TO = 1 defines the thermodynamic critical temperature 
Tc « 0.064593, it is larger than as expected. The 
value of the EA order parameter q at the glass transition 
is very close to 1, 1 — q = 0.00017116, so in this sense 
the phase transition is strongly first order but (similar 
to p-spin model), the energy and entropy do not change 
discontinuously at the transition. The mimerical solution 
shows that when the temperature decreases, the entropy 
of the glass state monotonically decreases and eventually 
becomes negative below T^ w 2.8 x lO^''. The explana- 
tion of such unphysical behavior is that IRSB ansatz, in 
fact, becomes unstable in this low temperature regime. 

Stability of the thermodynamical solution. In order to 
analyze stability of IRSB ansatz we expand the Eq.(|^) 
to the second order in fiuctuation of the order parameter 
5B and consider different families of fiuctuation matrices 
5B. This calculation is very similar to the analysis of the 
stability of paramagnetic solution and Parisi solution in 
SK model ||ll|,|l2| so we only sketch it here, for details see 
Appendix. We find that the most dangerous direction 
in the fiuctuation space corresponds to the "replicon" 
modes 1 11 12] that arc fiuctuations within diagonal blocks 
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oiSB satisfying the conditions {6B R)a,f3 = 0, 5Ba,a = 0. 
The eigenvalue corresponding to these modes is 

K = 2{1 - q)/ n + 2{l - qf - 2{r - q^) 

where r = J tanh^(z)P„i(z)(iz/ J Pm{z)dz. Numerical 
solution shows that Ai is positive at temperatures T > 
Tuns ~ 6.1 X 10~^ but changes sign at T^ns, thus IRSB 
solution is unstable at T < Tuns- 




FIG. 2. Main plot: Dependence of the configuration en- 
tropy on the temperature for the marginally stable solution. 
Insert: size of IRSB block, m, for thermodynamic (solid line) 
and marginally stable (dashed line) solutions. The value of 
T at at which m — 1 gives thermodynamical (dynamical) 
critical temperature. 



It is not clear however that the IRSB solution is a correct 
solution in the whole temperature range Tuns <T <Tg, 
on the contrary, it is quite likely that another solution 
is preferred by the system below some < Tg. We 
have only indirect arguments for this: first, it seems un- 
physical that Sconf decreases with temperature decrease, 
usually at lower temperatures additional states appear. 
Second, the obtained Sconf does not match the results of 
numerics if one believes that this solution remains correct 
at r < Tm- Finally note the analogy with higher tem- 
peratures: paramagnetic solution is always stable but is 
eventually replaced by IRSB solution. 

In conclusion we have justified the application of 
replica method to some systems without quenched dis- 
order and discussed situations in which cases it fails. We 
identify two dangers: correlations between metastable 
states close to the ground state and too large degeneracy 
of the ground state. We apply the formalism to the pe- 
riodic Ising spin model with orthogonal coupling matrix 
and find that it gives the same free energy as fiduciary 
Hamiltonian approach [Q . Further we show that it works 
in the intermediate temperature range but fails at low 
temperatures when metastable states become correlated. 
Two questions remain open: whether generalization of 
this method to continuous symmetry breaking would al- 
low one to study the models with correlated metastable 
states and what to do if the ground state of the model 
is highly degenerate as it is, e.g., in the case of unitary 
coupling matrix. 
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I. APPENDIX. STABILITY ANALYSIS. 



Here we provide the details of the stability analysis of IRSB solution. We start by expanding the Eq. 
second order in fluctuation of the order parameter SB: 

SF = -SBa.p Ma.Jl^jS SBj^s, 



to the 



where the coefficients a, b, c come from the expansion of Tr ^{B)/2 

a = 2q^- 2[Yrj/^i + q] [Yr]/{fi + 77m) - g], 



(17) 



(18) 



b^qifi 1 + 4(1 - q)] + Yri [(1 - q)rim - l]/[fj.^ + ri^xm], 



(19) 



c = 2(l-9) [2(l-q)+^i-% 



(20) 



with Y = 1 — g + qm. The function T is the four spin replica correlation function r^^^^-y^^ = (SaSpS-^Ss), defined with 
respect to the Hamiltonian (||). When at least two indexes, e.g. a and /3, are equal this function is 



^a,a,j,s — (1 ^ 9) Syj + qR^^s- 



(21) 



If a,/3,7,(5 are all different, the function V is not zero only if cither a,/?, 7, (5 belong to the same block, or two of 
indexes (e.g. a, [3) belong to one block and the other two (7, 5) to another one. In the first case Ta.p,^,s = r, where r 
is 



P„i{z) tanii^ {z)dz / / Pm{z)dz 



In the second case, since a, (3 and 7, S belong to different blocks T^^p^-y^s — ■ 
The eigenvalue equation for the matrix M is 

Ma.p,j,S SB^^s = A SBa,p. 

The eigenvalue equation can be simplified using that SB is a symmetric matrix. The resulting equation is 
[b - 2q{l - q)] {SB R + RSB)c,,p + {a- 2q^) {RSB R)c,,p - 2{Aq - 3r - 1) Sa,p SB^^^a. 
+ (3g^ - r)Ra-fS B-f^sR'f,sRs.p + 2(g - r) 2 Sa,p{RSB)a,a + Ra,pSBa.a + R/3,aSBp,p 
+ {q^ - r) ISa^pRp^-ySB-,^^ + 2Ro,,pSBa,p - 2R^^p{SB i?)„,<, - 2Ra^p{SB R)p^p 



(22) 



(23) 



-SaA^^B R)a,a - Ran^B^^^R^^p = [A + 2{1 - qf ~ c] SB. 



(24) 



This equation has a block-diagonal structure, therefore one can divide the fluctuation matrix SB into blocks SB 
of m X TO size and consider the fluctuations within each block independently. Moreover, all eigenvalue equations 
corresponding to diagonal (off-diagonal) blocks are equivalent. Therefore we are left with two cases: (I) fluctuations 
within an off-diagonal block and (II) fluctuations within a diagonal block. 

We begin our analysis of eigenvalues with the case (I) for which the eigenvalue equation is reduced to 



(6 - 2q{l - q)){SB£ + RSB)a.p + {a~ 2q^){£ SB£)a.,p = A'SB^^p, 



(25) 



where A' = A — c+2(l — g)'^ and £ is the m x m matrix with all elements equal 1. The first eigenvalue corresponds 
to SB satisfying {SB£)a.p = and it is 



A«=c-2(l-g)2 



(26) 



The second eigenvalue is 



5 



A^^^ =c + mb- 2Y{1 - q) 



(27) 



and it corresponds to 5B satisfying {5B£)a,i) 7^ 0, {£dB£)a,0 = 0. The last eigenvalue of type I corresponds to SB = £ 
and it is 



A^^^ =c + mG-2Y'^, 



where G = ma + 2b. 

The eigenvalues of type II satisfy the equation 



(28) 



£a,0 [2{q - r){6Ba,a + dB^,^) + 2{r - q^)[{SB£)a,a + {SBS)^,^] + (r - q^)TrSB 

-d^jj [2(4g - 3r - 1)(5B„.„ + 4(r - q){£ 6B)^^o. + (r - q^)[TY5B - {£ SB f 
(-2g(l -q) + b){5B£ + £ 5B)^,0 + {a-r + q^){£ 6B£)^^0 = [A' + 2(r - q^)]SBo,,0 (29) 

The first eigenvalue of this type corresponds to {SB£)a,/3 = and it is 

Af^ =c-2(l-g)2-2(r-g2). (30) 

All other eigenvalues of type II correspond to SB of the form SBa,0 = Xa + V/b £a,f3 + Va £a,i3, where Xa, Ua are such 
that either Xa = ?/„ = or aj^ = x, Ua = y- The corresponding eigenvalues are: 

^2^3 = c+[mb + 2E± ^J{bm + 2Ey - 16i;6]/2, 



hfl = c+[mG + F± ^{mG + - 4GF]/2, 

where E = Aq — ir — 1 + m{r — q) and F = 4[3^/9r7^]/(r?, to)/(1 — m)m. The numerical solution shows that for the 
thermodynamical solution all eigenvalues are positive for temperature higher than T^ns- At temperature less than T^ns 
the eigenvalue K\ corrc^sponding to the "replicon" mode becomes negative. The marginal ("dynamical") solution is 
defined by A^^-* = 0, all the other eigenvalues for this solution are strictly positive in the region where it exists. Thus 
the replicon mode is always the most relevant fluctuation. 
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